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AIRBORNE BALLISTIC CAMERA TRACKING SYSTEMS 


William L. Redish 
Langley Research Center 

SUMMARY 

On April 27, 1968, Reentry-F mission, a turbulence heating experiment (ballistic 
reentry), provided the opportunity for an operational test of airborne ballistic camera 
tracking systems and methods. The use of postflight data to develop reduction techniques 
for this type of tracking mode was also made possible. 

Two aircraft instrumented with ballistic cameras were flown as an independent 
tracking system during the mission. A well-established ground-based radar net was used 
for normal mission requirements and was assumed as the standard for the mission. The 
feasibility and overall accuracy of the airborne system were determined essentially by a 
direct comparison of data between the two systems. 

This paper describes the data requirements and the special data-reduction tech- 
niques developed. A comparison study between the ground-based radar system and the 
airborne system is also presented. 


INTRODUCTION 

Various orbital and reentry space tracking missions require tracking capabilities 
beyond the normal range of established land-based stations. These requirements indicate 
the importance of studying the feasibility of airborne or shipborne installations. 

Aircraft installations have an advantage over land and ship networks because the 
aircraft can rise above local inclement weather. In operations where tracking areas may 
be changed prior to launch, aircraft installations offer the flexibility to accommodate these 
changes. The studies contained herein could affect the future planning and consideration 
of other types of airborne tracking systems such as radar -bearing aircraft used in con- 
junction with real-time star tracking devices. 

Preliminary studies were undertaken at the National Aeronautics and Space Admin- 
istration (NASA) Langley Research Center to determine the requirements for aircraft 
tracking stations using ballistic cameras. The feasibility of various data-reduction tech- 
niques was examined and checked by computer simulations. These studies mdicated that 
mobile tracking stations presented special problems with regard to aircraft location, 



camera orientation relative to the star field, and a moving baseline between the aircraft. 
Such problems had to be solved to compute the time-height profile of a reentry target 
vehicle with sufficient accuracy for the test mission requirements of a ±150-meter 
uncertainty. 

The initial theoretical studies indicated that these problems could be minimized by 
the adaptation of special mathematical constraints. These constraints required specific 
measurements to be made aboard the airborne tracking stations. The required measure- 
ments onboard the aircraft are described in this paper. The data analysis procedures 
involving departures from standard fixed-camera methods are also described. The corre- 
sponding equations are discussed in the appendixes. These equations are outlined in gen- 
eral and are presented as guidelines for the future planning of similar missions. 

The Reentry-F mission (ref. 1) provided the opportunity for a full operational test 
of the data -reduction techniques and equipment development. This mission, powered by 
the Scout launch vehicle S-164, boosted a payload into a ballistic reentry trajectory for 
spectra analysis of a special heat shield. Ground-based radar data furnished a reference 
for feasibility comparison with other tracking methods. The results of this comparison 
are presented in the figures. 


SYMBOLS 


A-c> B c>C c 


^x> B x> ^x 


Ay,By,Cy 


B 

DE 

D jj°k 


D jk 


azimuth angle, station coordinate system 

direction numbers of plane containing target and both tracking stations 
distortion coefficients, x-axis 
distortion coefficients, y-axis 

magnitude of baseline between aircraft (See figs. 6(a) and 6(b).) 

difference between true zenith angle and refracted zenith angle 

distance between any camera -plate image j or image k and plate center 
measured within plane of plate (See fig. 11.) 

distance between two camera -plate images j and k measured within plane 
of plate (See fig. 11.) 
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D,j,j distance between target image and star image ] measured within plane of 

plate 

E elevation angle, station coordinate system 

f camera focal length (See fig. 11.) 

GHA Greenwich hour angle of Aries, west 

h geodetic height 

h geodetic height of primary aircraft (See figs. 6(a) and 6(b).) 

I - ' 

h geodetic height of secondary aircraft (See figs. 6(a) and 6(b).) 

L longitude, east 

jyjC 2 number of statistical equations for M stars taken two at a time 

N geodetic normal vector, primary aircraft (See figs. 6(a) and 6(b).) 

N geodetic normal vector, secondary aircraft, after mathematical adjustment 

b 

(See figs. 6(a) and 6(b).) 

Ng geodetic normal vector, secondary aircraft, based on navigator's position 

estimate (See figs. 6(a) and 6(b).) 

P primary aircraft position (See figs. 6(a) and 6(b).) 

Q computed target position (See figs. 6(a) and 6(b).) 

Q c fractional adjustment of surface refraction because of aircraft altitude 

R Earth radius at aircraft position (See fig. 10.) 

Rg Earth equatorial radius 

Rj star residual (radial displacement) 
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RM 


R Xi 


R yj 


g 

w T 


X a’ Y a’ Z a 


X b’\ ,Z b 

*r Y i 

X mj’ Y mj 


X P’ Y P 


V Y P C 
Vpj 

^k’ Y pk 


x t ,y t 


mean star residual 

star residual (x-axis component) 

star residual (y-axis component) 

distance from camera optical center to plate image 

mathematically adjusted estimate of secondary aircraft position (See fig. 6(b).) 

navigator’s estimate of secondary aircraft position (See fig. 6(a).) 

target line of sight from primary aircraft (See figs. 6(a) and 6(b).) 

target line of sight from secondary aircraft (See figs. 6(a) and 6(b).) 

epoch time of observations 

time of Greenwich hour angle 

Earth rotation rate 

Cartesian coordinates of primary aircraft, sidereal system 
Cartesian coordinates of secondary aircraft, sidereal system 
star residuals on x- and y-axes, respectively 
measured plate coordinates, uncorrected for distortion 
comparator axes 

coordinates of focal point on camera plate 
plate coordinates of star j, corrected for distortion 
plate coordinates of star k, corrected for distortion 
plate coordinates of target angle, corrected for distortion 
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Z zenith angle at station geodetic position, measured from geodetic normal 

vector 

Zjj scaled zenith angle to horizon at aircraft altitude (See fig. 10.) 

Z zenith angle of target vector from primary aircraft (See figs. 6(a) and 6(b).) 

tr 

Z zenith angle of target vector from secondary aircraft after mathematical 

adjustment (See fig. 6(b).) 

Z ' zenith angle of target vector from secondary aircraft before mathematical 

s 

adjustment (See fig. 6(a).) 

Z scale zenith angle for adaptation of surface refraction to corresponding high 

altitude effects 

Z t true zenith angle of star at station geodetic position 

Q'jjj angle subtended by plate images with vertex at camera optical center 

Oq altitude refraction coefficient 

AH difference in horizon zenith angle between surface and high altitude 

(approximation) 

e Earth structural ellipticity 

T right ascension (sidereal longitude) 

direction cosines of vector relative to the X,Y, Z axes, respectively, of a 
specified coordinate system 

apparent direction cosines of control star j in the sidereal system 
<p geodetic latitude 

Subscripts: 

a primary aircraft 
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app designates apparent (refracted wave front of light) 

b secondary aircraft 

c aircraft subscript 

j,k star image designations 

local station coordinate system with vertex of axes at station position, X-axis 

pointed toward polar axis of Earth, Z-axis oriented along geodetic normal 
at station position, and Y-axis forming right-handed coordinate system 

m designates measured plate coordinates 

sa apparent star sidereal angle 

sidereal designates sidereal coordinate system, with X-axis pointing to first point of 
Aries 

T designates target plate image 

x,y,z designates specific coordinate axes for rotation of angles from one coordinate 

system to another coordinate system 

AIRBORNE DATA ACQUISITION SYSTEM 

The test mission operation used two aircraft: NASA 432 and NASA 438. Special 
instrumentation was installed on each aircraft. The instrumentation requirements are 
listed in figure 1. 

Each aircraft contained a ballistic camera bolted to the aircraft viewport. Each 
camera had time -coded rotating shutters which produced an intermittent time -coded 
sequence of exposures of varying time durations. This system facilitated the time corre- 
lation between star -trail exposures and target-trail exposures on the individual camera 
plate. 

The airborne system required one aircraft position to be known accurately. This 
knowledge was to be acquired through a loran C geodetic positioning system with a loran 
receiver installed on the primary aircraft. A multichannel recorder (MCTS) was mounted 
on each aircraft. This recorder produced the data tapes for subsequent use in the data 
analysis. (Data channels for respective aircraft are listed in fig. 1.) 
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On the primary aircraft, a radar system was adapted for the dual purposes of 
measuring the aircraft altitude and the magnitude of the instantaneous baseline between 
aircraft. The measurements were accomplished in an alternating sequence. On the 
secondary aircraft, a radar altimeter was installed for the purpose of measuring aircraft 
altitude. A timing mechanism system on each aircraft recorded the event times for the 
recorder channels on the respective aircraft. There was no time synchronization between 
aircraft. 

Figures 1 and 2 illustrate the sequence of procedures from the data acquisition 
onboard the aircraft to the entry of data into the mam computer data-reduction program. 
The auxiliary data sequences are shown in figure 1. The camera-plate analysis sequence 
of procedures is outlined in figure 2. After the mission, the tapes from the MCTS were 
processed through the data-reduction sequence noted in figure 1. 

A separate computer program for the reduction of the loran C primary aircraft 
position data was obtained from the U.S. Coast Guard. Some minor modifications were 
made before the program was used with the mission data. 

AIRCRAFT FLIGHT MODE 

Ground -station ballistic cameras are fixed to the surface of the Earth, and the star 
images record only the motion of the Earth relative to the star field. For short periods 
of time, the star trails are normally separate and distinguishable. 

On the test mission, the ballistic cameras were bolted to the aircraft frame. There- 
fore, random aircraft motion could superimpose the separate camera exposures so that 
the individual time-coded exposures of the star trail could not be distinguished. During 
the short exposure time that occurs with reentry missions, it was sufficient to fly the air- 
craft in a gradual turning pattern. This pattern tended to spread these exposures across 
the camera plate so that each time exposure was normally separate and distinguishable. 

It is possible, however, that excessive aircraft roll, pitch, and yaw might interfere with 
this process by superimposing some of the chopped segments (overlapping one exposure 
segment on top of another). The merged exposures were ignored when reading the plate 
data. 

For diffused targets such as those in the barium ion cloud experiment (refs. 2 
and 3), this operational flight mode would not be practical because of the superimposition 
of image exposure of the cloud mass. This super imposition would cause erroneous den- 
sitometer readings on the cloud image, and thus prevent an accurate estimate of the mass 
center line. 
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TARGET AND STAR FIELD COVERAGE 


On the Reentry-F mission, the aircraft flew the turning pattern with the baseline 
essentially parallel to and south of the target track. Although the camera of each respec- 
tive aircraft covered the respective targets for a longer period of time than is included 
in the subsequent data analysis, the simultaneous time overlap was approximately 
2.35 seconds. The recorded time exposures for the two aircraft are presented in table I. 

As a result of aircraft motion, some exposures were superimposed and could not 
be separated when the camera plates were analyzed. Consequently, the effective time 
overlap was reduced to approximately 1 second during which time the plates could be 
read accurately. 

Plate 121224 (NASA 432) acquired 11 stars that were readable. Plate 121217 
acquired only five stars of sufficient brightness to be read and identified accurately. 

Plate 121216 (NASA 432) was exposed by a separate camera but was not used in the com- 
puter reduction process. The following section discusses plate 121216 as it illustrates 
the effects of aircraft motion on plate readings. 

CAMERA-PLATE DATA DESCRIPTION 

The effects of aircraft motion are noted in figure 3. A section of plate 121216 was 
selected and enlarged for this illustration. The camera shutters operated for a period of 
time prior to target acquisition. The small trace (star trail) illustrates the effects of 
aircraft motion on the resultant plate exposure during the total operating time of the cam- 
era. The eight target -trail exposures correspond to eight of the star -trail exposures 
(see fig. 3) and did not reflect appreciable aircraft motion during the short timespan of 
the target exposures in this specific enlargement. 

Figure 3 also illustrates an exposure problem. An acceptable exposure for the star 
trail (small chops) resulted in an overexposure of the target trail (large chops). This 
overexposure caused considerable difficulty when reading the magnified image on the 
stereo comparator. Reasonable accuracy was obtained by reading the break points at the 
corners of the rectangular target chops to determine the beginning or ending of the chop. 

A mean between the two long sides of the rectangle was read for the center line of the 
target trail. 

Plate 121216 (fig. 4) was selected to illustrate the effects of aircraft motion on the 
star plate readings. Interpolation between the fourth and fifth point on the graph could 
not be accomplished with the accuracy of the smooth curve segments on either side. 
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TIME CORRELATION OF DATA 


In the Reentry-F mission, the recorded data used were curve fitted and time inter- 
polated before entering the computer reduction system. This procedure time -synchronized 
the data prior to the machine reduction process. The data were then entered in the com- 
puter program as if the time -code cycles of the camera shutters had been synchronized 
between the two aircraft during the mission. This method required extensive manual 
operation and greatly increased the operational time required for the reduction. Because 
of the experimental nature of the problem, this procedure was chosen for the following 
reasons: 

(1) The overall complexity of the computer program was reduced so that all primary 
reduction procedures could be debugged and optimized. 

(2) All input data were monitored manually for unexpected problems in the new proto- 
type system. 

(3) Machine curve fit error contributions were eliminated in the data output residuals 
during the development process. 

DATA ANALYSIS TECHNIQUES 

A general outline of the mathematical techniques involved in the computer reduction 
program is given in appendixes A to H. The data -reduction sequence adapted for the air- 
borne tracking system is illustrated in the flow diagram of figure 5. In figure 5, loop 2 
and loop 3 control mathematical iteration cycles. Each camera plate was processed sepa- 
rately. Loop 1 represents the return for the second camera plate. 

Astronomical star data for a specific epoch are compiled in star catalogs. The com- 
piled star data must be updated to the tracking date and time. This procedure is briefly 
discussed in appendix A. 

Since the camera plates were located within the atmosphere, the star images on the 
plates have been refracted. The astronomical direction cosines (true direction angles) 
must be transformed to refracted (apparent) direction cosines to be consistent with the 
observed plate data. The process used is described in appendix B. 

Initially, the focal length and camera-plate center were not known accurately enough 
because of the unknown factors of lens and aircraft viewport distortion characteristics. 
These parameters must be solved for by a least-square matrix iteration process during 
the camera orientation. (See appendixes C and D.) During each iteration cycle (loop 3), 
the data editing and weighting procedures were also calculated. (See appendix E.) 
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After the residuals were minimized for the respective camera plate, the target 
angle for that plate was computed by the equations of appendix F. The entire process was 
then repeated using loop 1 for the second camera plate. 

After the target angle for each plate was computed, the secondary aircraft was 
positioned with respect to the primary aircraft by the mathematical process described 
in appendix G. This geometric procedure is illustrated in figures 6(a) and 6(b) and is 
described as follows : 

The position of the primary aircraft P was obtained by a loran C system. This 
position was assumed to be correct in the data- reduction process. An initial estimate 
of the secondary aircraft position s' was taken from the onboard navigator's data 
(fig. 6(a)). The adjusted secondary aircraft position s (fig. 6(b)) was subsequently deter- 
mined by a mathematical process using the following data. 

The magnitude of the baseline between the two stations B was measured by a 

radar -ranging system mounted on the aircraft. The respective geodetic heights h 

and h g were determined by aircraft radar altimeters. The orientation of the target 

lines of sight T^ and T g was determined by the ballistic cameras mounted on the 

viewports of each aircraft. This orientation used the star background for control points; 

the subsequent orientation angles were measured in the sidereal (star) coordinate system. 

The vectors N p , N g , and N g represent geodetic normals to the theoretical ellipsoid 

of the Earth at the respective geodetic latitudes and longitudes of the respective stations. 

The angles Z. Z , and Z' represent the zenith angles computed to estimate refrac- 
p s s 

tion corrections of star and target data. 

In the reduction process, the initial estimated position of the secondary station s' 

was improved by forcing the target lines of sight T and T to intersect at Q while 

s p 

the measured range B and the altitudes h and h continued to be maintained. It 

p s 

should be noted that the target lines of sight T and T are rigidly determined by 

p fa 

camera -plate reduction and that these lines of sight may be freely translated geometri- 
cally without changing their orientation. Therefore, when the initial secondary aircraft 
position s' is mathematically adjusted to a more correct position s the line T g is 
translated without changing its orientation. The two lines T and the displacement are 
parallel because they represent the same orientation in the specified coordinate system. 
This process is repeated at each data point time epoch. After the secondary aircraft 
position was adjusted, the secondary aircraft camera plate was reprocessed through loop 2 
of figure 5. The processing of the primary aircraft was negated after the first cycle of 
computations m loop 2. The final triangulation for target height and position was computed 
in the sequence as noted in figure 5. The selected process is described in appendix H. 
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DATA ANALYSIS RESULTS 


Camera-Plate Processing 

Ballistic camera plate 121224 produced 11 positively identified stars that were 
readable. The readable stars were adequately displaced over a sufficient plate area to 
facilitate distortion analysis. The iteration procedure for computing the distortion coef- 
ficient values reduced the mean star residual errors from 17.3 pm to 4.8 pm. This final 
value was close to the 2-pm value normally expected with ground -station fixed-camera 
systems. 

Ballistic camera plate 121217 produced only five stars of sufficient magnitude to 
be read and identified. The small number of stars prevented distortion analysis for this 
plate. As a result, a statistical method based on successive approximations of the plate 
star residual errors was adapted to compensate partially for the camera and aircraft 
viewport distortion errors. The mean star residuals were improved from an initial 
56. 1 pm to a final 29. 1 pm. The grouping of the stars on one side of this plate probably 
caused the higher initial mean star residual error as compared with the initial mean 
error of plate 121224. 


Triangulation Processing 

Table n lists the results of the triangulation for the succeeding time epochs com- 
puted. The time -height profile graph for the payload target is shown in figure 7. The 
solid line represents the statistically smoothed ground radar track and was used as a 
basis of comparison. The circles represent airborne camera computed heights with no 
refraction corrections (minimum). The triangles represent the computed heights with 
surface refraction corrections (maximum). The final target computed heights with the 
assumed high altitude refraction corrections are represented by shaded diamonds and 
occur between the two computations shown. The height variation caused by the refraction 
uncertainty appears to be less than the 150 -meter scale indicated on the graph. 

Figures 8 and 9 illustrate the target position (latitude and longitude) relative to the 
computed position. The loran C equipment was operating at its maximum effective dis- 
tance for reasonable accuracy. This distance probably resulted in the bias which offset 
the entire tracking net computed position as indicated on the graph. However, since the 
ranging system between aircraft maintained a continuous and accurate estimate of rela- 
tive baseline distances, the time -height profile was not significantly affected. 
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RECOMMENDATIONS FOR DATA IMPROVEMENTS 


Although the feasibility of the airborne system time -height profile accuracy has been 
proven under marginal conditions, these conditions are not desirable for normal range 
operations. The marginal conditions are discussed in this section. 

The aircraft viewport is a curved surface when viewed with the accuracy of a ballis- 
tic camera. During flight, the aircraft undergoes vibration and structural bending which 
cause variable window distortion and subsequent optical effects upon the entering wave 
front of light. A larger number of control stars well dispersed over the entire plate are 
needed for better resolution of these distortion effects. A minimum of 15 control stars 
is recommended to facilitate the editing and rejection of erroneous control-star data. 

Three recommendations for acquiring a larger number of control stars are: 

(1) Incorporate image intensification techniques to enhance the readability of dim 
stars without overexposing the target trail. 

(2) Plan launch window constraints based on observable star field densities. 

(3) Acquire camera-plate exposures over high star density areas during the aircraft 
approach to the mission tracking area. 

The last suggestion would allow a reasonable distortion and refraction analysis in the event 
that the subsequent mission track covered a darker portion of the sky. 

The use of time synchronization between camera shutters would tend to minimize 

the data -reduction process. This procedure would eliminate the necessity for adapting 

/ 

preliminary manual or machine data -interpolation procedures prior to the main reduction 
process. 


CONCLUDING REMARKS 

An airborne ballistic camera tracking system was tested during a ballistic reentry 
space mission and was compared with a well-established ground-based tracking radar net- 
work. The experimental airborne system computed a time-height profile that compared 
favorably with the results from a ground-based radar system. 

Airborne ballistic camera -plate data with an adequate number of recorded stars well 
dispersed across the entire plate may be reduced to a 5-fj.m plate accuracy. The final 
time -height target profile between the airborne system and the established ground system 
differed in altitude by less than 150 meters despite some marginal data. 
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Since this profile was within the limits initially set for the project, airborne tracking 
networks appear feasible. Some areas of necessary improvement to msure the accuracy 
of future missions are presented. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
May 4, 1976 
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APPENDIX A 


STAR DATA PROCESSING 

The star data of the identified plate stars were obtained from a standard star catalog 
for the 1950 epoch (ref. 4). This information was updated to the mission time in the com- 
puter program. 

Firsts and second-order polynomials were found to reproduce the star precession 
tables to the eighth significant figure. Since this calculation was adequate for a century 
of updating, these coefficients were built into the computer program. The nine preces- 
sional parameters were computed from these coefficients and were used in a rotational 
matrix to correct the star data. The star proper motion was also corrected to the cur- 
rent mission epoch. 
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APPENDIX B 


REFRACTION CORRECTIONS 


The corrected star angles (from appendix A) are true angles (within astronomical 
accuracy). Since the camera plates were exposed in the atmosphere, the plates actually 
record the refracted wave fronts of light. The star angles must be corrected for refrac- 
tion before they can be used as orientation control points for a mission. 

Since refraction is a local phenomenon in the vicinity of the tracking station, the 
Earth-centered sidereal star direction cosines must be rotated to a local coordinate sys- 
tem with the tracking station as the center of the system. The local system was defined 
with the x-axis pointing toward the north polar axis. The z-axis pointed along the local 
vertical, and the y-axis pointed west to form a right-handed coordinate system. The 
rotation is shown below in matrix form: 


X local 



^x 

"x~ 


fx 1 

sidereal 

^ local 

= 

x y 


v y 

• 

^sidereal 

"local 


X z 

w 

^z 

"z 


"sidereal 


(Bl) 


The nine directional cosines X x to i> z are obtained from the best estimate of the 
station latitude and longitude. If the station position is adjusted during the reduction pro- 
cess, the apparent star angles should be recomputed before any further analysis is 
undertaken. 


Geodetic longitude is changed to right ascension (sidereal longitude) by the relation- 
ship where 


r= L + GHA + W £ (t - t ) 

L station geodetic longitude, east 


(B2) 


T right ascension of geodetic normal 

GHA Greenwich hour angle of Aries (from Air Almanac) measured west of 

Greenwich meridian 


Wg Earth rotation rate 
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APPENDIX B 


t time of current epoch 


t 


g 


time of Greenwich hour angle, west 


The nine rotational matrix elements are then computed from the foregoing declina- 
tion and the right ascension of the geodetic normal as follows. 


X = -sin (p cos r 
P x = -sin cp sin r 
v = cos ( p 
Ay = sin r 
Py = -cos r 



A z = cos (p cos T 
P z = cos (p sin r 
i> 2 = sm (p 


(B3) 
(B4) 
(B5) 
(B6) 
(B7) 
(B8) 
(B9) 
(BIO) 
(B 1 1) 


The parameter <p is the declination of the geodetic normal (in this case, station geodetic 
latitude). The rotation into the local system is thus completed. 

The true and the apparent direction cosines of the star vector of orientation may be 
expressed in terms of station azimuth A and elevation E. 'The azimuth is unchanged 
by refraction. The elevation angle (derived from the cofunction of the zenith angle mea- 
sured from the geodetic normal) has both true and apparent relationships. These are: 


"local = cos E cos A 

(B12) 

"local - cos E sin A 

(B 13) 

local = sln E 

(B 14) 

"app = cos E app cos A 

(B15) 
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APPENDIX B 


= cos E„__ sin A 
^app app 


v - sin E 
app app 


(B 16) 
(B17) 


where the subscript app denotes the apparent (refracted) angle. 

The true elevation E and azimuth angles A were solved for in equations (B12) 
to (B14). The apparent elevation angle E ap p must be computed from E with some 
system of refraction corrections. The system used is described in the following 
paragraph. 

The elevation angle E may be obtained from equation (B14) after the matrix rota- 
tion (eq. (Bl)) is completed. The angle becomes 


E = arc sin ^ local 


(B18) 


The apparent direction cosine relationships are obtained by combining equations (B12) and 
(B15) and by combining equations (B13) and (B16) by the azimuth terms. Equation (B17) 
is unchanged: 


X app 


cos E 
cos 


app 

E 


* local 


(B19) 


cos E 


M 


app 


app 


cos E M local 


(B20) 


v 

app 


sin E 


app 


(B21) 


An interim system based on a modification of reference 5 was adapted for the test 
mission. Additional information and techniques concerning high-altitude refraction may 
be obtained from references 6 to 8. 


The surface refraction correction DE 
zenith angle Z and the true zenith angle 

“rr 


DE = Z - Z 

sc app 


(ref. 5) is correlated with the apparent 

as shown: 
sc 


(B22) 


For higher altitudes, the value of DE decreases. This decreased value was modified 
by assuming an exponential scale gradient for a coefficient of DE. The scale gradient 
adapted is ' 
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APPENDIX B 


Q 


c 


-o-oh 

e 


(B23) 


where Q c represents the fractional change of DE because of altitude h. When equa- 
tion (B23) was adapted as a coefficient in equation (B22) and was rearranged, the first 
modification became 


Z a pp = Z sc - Qc< DE > < B24 > 

The exponential coefficient a Q was experimentally determined by successive trial 
and error computations using plate 121224 with 11 stars. The accepted value was based 
on minimization of the mean star residuals on this plate. 

A second modification was adapted and contributed to the minimization of the plate 
star residuals. The parameter DE (described in ref. 5) is based on a true zenith hori- 
zon angle of 90° (tt/2 radians). The parameter Z represents the true zenith angle Z f 
of the star. At the aircraft altitudes flown during the test mission, the zenith angle of the 
true horizon is slightly greater than 90° (illustrated in fig. 10). The parameter Z then 
becomes a function of the true zenith angle Z^. A proportional scale modification based 
on figure 10 was assumed. This modification is described in the next paragraph. 

A local Earth radius R was computed for the specific aircraft position. From the 
geometry of figure 10, the maximum increase in zenith angle AH due to station altitude 
was approximated by the following geometric relationship: 

AH « arc cos ( =3 — ) (B25) 

\R + hj 

where h is the aircraft geodetic height. The increased horizon zenith angle Z^ at the 
station altitude h is 

Ztt = ^ + arc cos ( _ R ) (B26) 

n Z \R + h / 


The scale ratio is 



(B27) 


18 



APPENDIX B 


When equation (B26) is substituted into equation (B27) by Z H and the result is substituted 
into equation (B24) by Z , the apparent zenith angle becomes 

SC 


IW 


app 


5 + 


arc cos 


R 


R + 


h/J 


- Q r (DE) 

\1 c 


(B28) 


If the zenith angles Z f and Z„ of equation (B28) are expressed m terms of their 

L ^-PP 

corresponding elevation angles E and E„ : 

app 




(B29) 


and 


Z = - - E 
app 2 


app 


the relationship becomes 


E app 


l (AH + E) 

~ + Qc< DE > 

2 + AH 


(B30) 


(B31) 


The parameters Qg and AH are the relationships expressed by equations (B23) 
and (B25), respectively. Equation (B31) may be substituted into equations (B19) to (B21) 
to compute the apparent direction cosines A app , M app > and y . 

These apparent star angles may then be rotated back to the sidereal system by an 
inverse matrix rotation as shown: 


X sa 


X x *y X z 


X app 

(B32) 

CG 

=1 

= 

^x ^y ^z 


^app 

(B33) 

u sa 


"x ^y 


v 

app 

(B34) 


The subscript sa designates the apparent star sidereal angle exposed on the ballistic 
camera plates. These final direction cosines were then used in the following procedures. 
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CAMERA ORIENTATION EQUATIONS 


In previous operations, fixed-camera systems have usually been oriented to a geo- 
detic system involving a rigidly surveyed local vertical. The interior camera geometry 
and the star field were related to the local vertical by initially using the best estimate of 
the station position and solving for the following parameters in a least-square matrix 
iteration. 

(1) Focal length of camera; 

(2) Plate center (intersection of the camera focal line with the plane of the exposed 
plate or film) in the plate coordinate system; 

(3) Camera tilt angle about the focal line; 

(4) Right ascension and declination (or azimuth and elevation) of the focal line vector. 

Aircraft motion and the complexity of stabilized platforms indicated that a new 
method which excluded the local vertical was needed. It was noted that the star field 
(projected upon the camera plate) constituted a unique set of geometric interrelationships 
when projected through the camera geometry. Relative star displacements on the exposed 
camera plates (or film) represented relative angular displacements of the stars projected 
onto a celestial sphere. Although the exposure was made near the surface of the Earth, 
these angles are essentially equal to those measured from the center of the Earth because 
of the absence of measurable star parallax. Because of this similarity, the measured dis- 
placement angle between two specific stars is essentially constant from any pomt in the 
solar system and, therefore, independent of any planetary geodetic coordinate system. 

This criterion indicated that a basic cosine equation was adequate to describe the relative 
displacement angles of points on the plate. 

The basic equation also excluded the necessity of solving for the pointing angle of the 

focal line vector and the camera tilt angle. Three variables were thus eliminated from 

,/ 

the least-square matrix solutions, and a rapid convergence for the three remaining 
unknowns was insured. This equation and its related forms are described m the following 
paragraphs. 

The following general cosine equation is based on the inverted coordinates of the 
camera geometry of figure 11. This equation is used statistically to combine star 
pairs j and k over all available stars on the plate to minimize the residual errors of 
star images. The basic relationship is 




2r j r k COSO! jk 


(Cl) 
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where 



plate distance between plate star images 


r distance from camera optical center to a respective subscripted star image j 

or k on plate 


The angle with the vertex at the camera optical center is obtained by the dot 

product of any two subtending vectors with orientation X, /i, and u through correspond- 
ing images on the plate. The equation is given as follows: 


cos “jk = Vk + k + u )\ 


(C2) 


The parameters r, and r k represent the vector magnitude from the camera 
optical center to the respective plate images. These parameters relate to the focal 
length f and force the focal line vector to remain essentially perpendicular to the plate 
by the following relationships: 


f =f2 + D l 2 

(C3) 

k= f2 + D £ 

(C4) 


where 


D 


camera focal length 

distance from focal point on plate to respective subscripted star image on 
plate 


The distances (in the plane of the plate, between plate images) are computed from 
the following parameters: 

Xpj,Yp. arbitrary plate coordinates of star image j corrected for distortion 

Xpk ? Ypk arbitrary plate coordinates of star image k corrected for distortion 

Xpc’Ypc arbitrary plate coordinates of focal point on plate 
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The arbitrary plate coordinate system essentially constitutes the axes of the stereo- 
comparator plate or film reader. Since the basic cosine equation expresses the displace- 
ment differences, the actual location of the comparator axes Xp and is immaterial. 

In figure 11, the distances (normally corrected for distortion) between plate images 
are expressed as 


D i - * ( Y pi - V) 2 

(C5) 

^ = (V-V ) 2 + ( Y Pk-V ) 2 

(C6) 

i) 2 + ( Y pk - Y pj) 2 

(C7) 

When equations (C5) and (C6) are substituted into 
the following equations are obtained: 

equations (C3) and (C4), respectively, 

r f = f2 + (\i ' V ) 2 + ( Y Pi ' Y P c) 2 

(C8) 

4 = f 2 * (Xpk - Xpc) 2 + ( Y pk - vf 

(C9) 


The final equation results when the relationships of equations (C2), (C7), (C8), and (C9) 
are substituted into equation (Cl). This equation becomes 


(*pk - Xpj ) 2 + (v - Y pj ) 2 = 2f2 + (*W - v ) 2 + ^ - v ) 2 + (v - v*' 


+ (Y pk-V 


:) 2 - 2 /( j 


r + D 2 )(f 2 + Dg/XjXfc. + PjP k + Yk) 


(CIO) 


where Dj and are obtained from equations (C5) and (C6), respectively. 

The parameters 

after statistical correction for distortion. The relationship between the corrected coordi- 
nates and the actual plate measurements is discussed in appendix D. 

The combinations of star pairs produce a maximum statistical use of the star data 
available. For M stars taken two at a time, the number of statistical equations available 
are 


, Ypj, Xpp., and Y pk represent the plate image coordinates 


r = Ml 

M 2 2(M - 2) ! 


(Cll) 
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Therefore, for the five -star plate 121217 there are 10 statistical equations. The number 
of equations probably explains the statistical improvement of the accuracy of this plate as 
noted in the section "Data Analysis Results." 
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DISTORTION COEFFICIENT FORMS 


Numerical analyses indicated that for the Reentry-F mission, the distortion charac- 
teristics of the camera lens and aircraft viewport approximated a sine wave function on 
each axis. The subsequent relationships between the corrected plate image coordinates 
and the measured image coordinates follow: 


where 


Xp, » Ax sin |B x (x m] - Xp C ) - B x c x ~i + x mj 


-PJ 


= A v sin 1 

B (Y • - Y \ 

- B C 

y 

y\ mj pc) 

y y 


+ Y mj 


(Dl) 

(D2) 


V Y Pi 

Vmj 


A R C 

y’ r y 


j 


plate star image coordinates statistically corrected for distortion 

measured plate star image coordinates 

distortion coefficients, x-axis 

distortion coefficients, y-axis 

star images j or k 


If equations (Dl) and (D2) are substituted into equation (CIO), a final equation is 
obtained for each star pair processed. The unknowns f, Xp C , Yp C , B x , C x , Ay, 

By, and C y may then be solved over all available stars by a matrix iteration least-square 
process. The resulting plate residuals may then be computed directly from equations (Dl) 
and (D2) as shown: 

Residual (Xj) = X pj - X mj (D3) 

Residual (Yj) = Y pj - Y mj (D4) 

The parameters Xpj and Ypj represent the computed plate coordinates for star j, 
statistically corrected for distortion. The statistical mean residual over all stars was used 
as a criterion for determining the feasibility of the plate data-reduction processes. 
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DATA EDITING AND WEIGHTING PROCESS 


In the test mission, the star field for a specific plate presented the only control 
available for computing the related camera orientation. Erroneous control star data 
could conceivably have a significant effect on the final accuracy of the system. A rela- 
tive weighting process was therefore adapted for the purpose of down weighting incon- 
sistent star data within the field of view of each camera plate. 

The weighting process was based on the final residuals of equations (D3) and (D4) 
after the best estimates of the distortion coefficients were obtained. The computed 
weights for each available control star were used in the subsequent iterations (to loop 3, 
fig. 5) for statistically improved values of focal length, plate center location, and distor- 
tion coefficients. 

A Gaussian distribution form was adapted for this process. The radial displacement 
residual of each star was computed as shown: 

R i-K* R y, (E1 > 

where 


R^ residual Xj 

Ryj residual Y. 

as noted in equations (D3) and (D4). A statistical mean residual (RM) over all available 
stars on the plate was computed. The resultant weight was computed for each control 
star j by the following relationship: 


Weight (j) = 



(E2) 


These weights were subsequently changed for each iteration cycle of loop 3 (fig. 5). When 
the plate residuals were minimized through the overall iteration process, the final data 
were used for the computation of the target angle. 
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TARGET ANGLE REDUCTION 


After the best statistical estimates of the focal length, plate center, and distortion 
coefficients have been computed, any direction cosine of a target image on the plate may 
be computed by rearranging equations (Cl) and (C2) as shown: 


where 


and 


Xj(X T ) + Uj(n x ) + •'j(.'T) - 


2 2 
r T + r i 


D 


T] 


2r T r j 


(FI) 



x t ,y t 


Vpi 

V Y pc 
f 


plate target coordinates, corrected for distortion 
plate star coordinates, corrected for distortion 
coordinates of focal point on plate 
focal length of camera 


Arp apparent direction cosines of target image in sidereal system 

A.,u apparent direction cosines of control star in sidereal system 

x^J rJ rj 

This process produces one equation for each control star j. The target angles are 
unknown and may be solved for over all usable stars on the plate in a least-square sense. 
The resultant direction cosines A,p /irp and are then corrected for refraction and 
subsequently used for adjusting the secondary aircraft position and computing the final 
triangulation position as described in appendixes G and H. 
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RELATIVE POSITIONING WITHIN THE AIRCRAFT NET 

The relative positions of the aircraft within the net were determined in the actual 
data- reduction process. The process is illustrated in figures 6(a) and 6(b), The air- 
craft positions were measured in the geodetic coordinate system based on the Greenwich 
meridian. The orientation of the target lines of sight T and T were reduced in the 
sidereal coordinate system. Coordinate rotations were made from the geodetic to the 
sidereal system for compatibility of data. 

The target line of sight T was translated mathematically until it intersected with 

s 

line of sight Tp at point Q. Point Q is essentially correct when the measured dis- 
tance B and the measured height h are maintained during the translation process. 

The position s then represents the corrected relative position of the secondary aircraft, 
and the subsequent triangulation was based on this position. 

It should be noted that this method precludes residual analyses from the resultant 
triangulation. In this case, triangulation residuals merely indicate the mathematical 
accuracy of the entire boundary solution process. Consequently, a mission comparison 
with a well-established ground-based trackmg net was necessary to establish an estimate 
of the overall accuracy of the system. 

The mathematical adjustment is accomplished by the combination of the three non- 
linear constraints into a 2 by 2 determinant. The three constraints are described as 
follows : 

The first constraint concerns the primary aircraft position with subscript a and 
the measured range B to the secondary aircraft position with subscript b. This process 
is mathematically described by a sphere of radius B with the center at the primary geo- 
detic position: 



Theoretically, the two stations and the target point lie in the same plane when the 
secondary position is adjusted and the target lines of sight are forced to intersect at any 
time epoch condition. This instantaneous epoch plane (from ref. 9) can be expressed in 
linear mathematical form by 

(*b - X a) A c + ( y b - Y a) B c + - Z a) C c = 0 < G: 
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where 

A c = Va " "a^b 


(G3) 


B c = Vb 


XvV 

b a 


(G4) 


C c = Va “ Vb 


(G5) 


and X, ji, and v represent the target direction cosines of the respective aircraft. The 
respective direction cosines X, p, and v are the true (unrefracted) target directional 
cosines. 

The third constraint is expressed in Cartesian parametric form. These relationships 
contain the geodetic relationships required and are described as follows: 


Re 

Ke = - 

(G6) 

\] 1 - s in2 (j) 


x c = ^h c + Kg) cos (p cos r 

(G7) 

Y c = ^h c + K E ) cos (p sin Y 

(G8) 

Z c = [ h c + ( 1 - £2 ) K E] sin< '’ 

(G9) 


where 

h geodetic height of station above ellipsoid of Earth 

e structural ellipticity of ellipsoid of Earth 

Rg structural equatorial radius of ellipsoid of Earth 

T right ascension of station position 

c aircraft subscript a or b 

The right ascension r is obtained from the relationship in equation (B2). If equa- 
tions (G7), (G8), and (G9) are substituted into the constraint equations (Gl) and (G2) by the 
appropriately subscripted X, Y, and Z, two equations in two unknowns are available. 
The two unknowns are the latitude <f> and the longitude L of the secondary station. 
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When the two resultant nonlinear equations were linearized by Taylor's Series and formed 
into a 2 by 2 determinant, they were found to be highly convergent with five to seven itera- 
tions. The initial trial values for the iteration were not critical. 

This adjustment changed the orientation of the station geodetic normal from N’ of 

s 

figure 6(a) to N of figure 6(b). This process changes the zenith angle Z' to zenith 
s s 

angle Z ; subsequently, the resultant refraction correction for the camera plates of the 
s 

secondary station is changed. The change in refraction was between 1 and 2 percent of 
its prior value. Although this change is small, the entire plate reduction process was 
recomputed to obtain the maximum possible accuracy from the system. This adjustment 
completed the final data required for the resultant triangulation for target position and 
output parameters as described m appendix H. 
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TRIANGULATION 

After the secondary aircraft position was adjusted, the final target height and posi- 
tion were computed. The process selected was based on the geometric iteration illus- 
trated in figure 12. 

The iteration procedure is initially started from the position of the primary sta- 
tion. A line was projected from point p to point A and perpendicular to the fixed line of 
sight T . Another perpendicular was dropped from A to K. The iteration was continued 
until the minimum distance between T^ and T g was obtained. This minimum distance 
determined the convergence criterion. The final target point is assumed half-way between 
the two lines of sight along this minimum perpendicular distance if any error of intersec- 
tion remains. In the mobile tracking net, when the lines of sight are forced to intersect, 
no triangulation error can be determined. In fixed-station tracking nets, this minimum 
distance would represent a triangulation residual since the lines of sight are not forced 
to intersect by the adjustment procedure used in airborne systems. 
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TABLE I.- TIMING OVERLAP BETWEEN CAMERAS 


Station 

438 


Payload 


Plate 

121217 


Station 

432 


Plate 

121224 


Time, s 


Time, s 


448.91 

449.43 

449.96 
450.49 
451.02 
451.52 
452.08 
452.63 
453.13 
454.18 
454.70 

455.97 
456.29 


457.34 

457.88 

458.40 

458.93 


456.58 

457.10 


457.73 


.K rt 


o 

a> 


458.27 W 


458.79 


> 

o 


459.30 

460.27 


a Time is seconds from launch; read at end of each exposure. 


TABLE, II.- COMPUTED TARGET DATA 


Range time, s 

Height, km 

Latitude north 

Longitude west 

457.8 

20.283 

33.2402 

-62.9741 

458.0 

19.871 

33.2353 

-62.9635 

458.2 

19.437 

33.2312 

-62.9528 

458.4 

19.010 

33.2253 

-62.9422 

458.6 

18.591 

33.2197 

-62.9313 





r 


, 

Primary aircraft (A) | 



CAS Ballistic camera system (shutter timing) 
LRC Loran C time difference system 
MCTS Multichannel tape system 
RAS Radar altimeter system 

RRS Radar ranging system 

TMS Timing mechanism system 

Figure 1.- Auxiliary data processing. 
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Figure 2. - 



-plate data processing. 





















Figure 3.- Enlargement section, plate 121216, aircraft 432. 
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Figure 4.- Plate 121216, star No. 11593, XY -plate reading (p. m). 



Loop 1 



Figure 5.- Data -reduction sequence. 
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(b) Adjusted secondary aircraft position. 
Figure 6.- Secondary aircraft positions. 



Height, km 
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Time, sec from launch 


Figure 7. - Time -height profile. 
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Latitude, deg north 


33.32 




O 


33.30 


33.28 


33.26 


33.24 


33.22 


33.20 


33.18 




457.2 457.4 457.6 457.8 458.0 458.2 458.4 458.6 458.8 459.0 


Time, sec from launch 
Figure 8.- Time -latitude profile. 
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